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Abstract 

The inclusion of super symmetric threshold corrections to the 6-quark mass has dra- 
matic consequences in scenarios with large CP-mixing effects in the Higgs sector. In 
particular, when the phase of the combination M g u is ~ 180°±30°, the lightest sbottom 
squark becomes tachyonic and, possibly, the 6-quark Yukawa coupling nonperturbative 
for values of tan (5 ranging from intermediate up to large or very large, depending on 
the size of arg(Affj,), arg(Af ) fi), and the details of the spectrum. For these scenarios, 
when allowed, as well as scenarios with different values of arg(Mgu), the cross sections 
for the production of the three neutral Higgs bosons through 6-quark fusion have inter- 
esting dependences on arg(A t u) and arg(Abfi), and the deviations induced by the nib 
corrections are rather large. In general, such production channel cannot be neglected 
with respect to the production through gluon fusion. For large CP-mixing effects, the 
lightest neutral Higgs boson can be mainly CP odd and the 6-quark fusion becomes its 
main production mechanism. Searches at the Tevatron and the LHC can easily detect 
such a Higgs boson, or constrain the CP-violating scenarios that allow it. 



It is well known that threshold corrections to the fe-quark mass due to the virtual ex- 
change of supersymmetric particles can be large and do not decouple in the limit of heavy 
superpartners [T]. They turn out to be quite substantial |2], also for moderate values of tan j3, 
in scenarios that maximize the CP-mixing effects induced in the Higgs sector at the quantum 
level by phases in the supersymmetric and soft supersymmetry-breaking mass parameters 
(such scenarios were first studied in Refs. [Hill]). 

These corrections affect the vertex b-b-Hf, which is responsible for one of the production 
mechanisms of neutral Higgs bosons at hadron colliders. If no CP-mixing effects are present 
in the Higgs sector, Hf is either one of the two CP-even states h and H, or the CP-odd 
state A [5]. In this case, the mechanism of 6-quark fusion has relevance, for example, for the 
production of A 6j (and in the decoupling limit, also in the case of the production of H). 
Indeed, a factor tan (3 in the b-b-A coupling can give a production rate comparable to, or 
even dominant over, the rate of production through gluon fusion. On the contrary, lacking 
this tan f3 factor, the 6-quark fusion plays a minor role in the production of the light CP-even 
Higgs boson with respect to the gluon-fusion, which, although loop-induced, is advantaged 
by the large yield of gluons in the proton. In scenarios with a CP- violating Higgs sector, the 
three mass eigenstates are mixed states with both CP-even and CP-odd components *. 
Thus, it is conceivable that, if the CP- violating mixing effects are large, the production 
through 6-quark fusion may be relevant for all mass eigenstates of neutral Higgs bosons. 

The aim of this paper is to investigate this issue, while consistently including the super- 
symmetric threshold corrections to the 6-quark mass. We find that these corrections can 
have a strong impact if phases are present in supersymmetric and supersymmetry-breaking 
parameters, in that they can exclude intermediate/large values of tan/3, leaving nevertheless 
allowed the very large ones. Predictably, when these corrections and the CP-mixing effects 
in the Higgs sector are large, the production cross sections through 6-quark fusion for all the 
neutral Higgs bosons are affected, with that for the lightest one, possibly, dramatically en- 
hanced. The effect of the supersymmetric threshold corrections to m;, on these cross sections 
may remain rather large also when some phases in these scenarios acquire trivial values, i.e. 
0° or 180°. 

We illustrate our findings for scenarios that have been dubbed CPX scenarios JUj; and 
that tend to maximize the CP-mixing effects in the Higgs sector *. These scenarios are 
identified by the spectrum: 

\A t \ = \A b \ = 2c A Msusy , M = 4c M M SUS y , rn^fi, = M SVSY , \M~\ = 1 TeV , (1) 

where A t and Ab are the complex trilinear soft terms for the third generation squarks, /i is 
the complex supersymmetric Higgs(ino) mass parameter; ^q 3 , m ij 3 an d m D 3 are t ne rea l 

*Phenomenological consequences of these mixings can be found in Refs. [Z| 

Hn general, large CP-violating phases in supersymmetric models are indirectly forbidden by the non- 
observation of electron and neutron electric dipole moments. These constraints, however, can be evaded 
by cancellations between the one- and higher-loop contributions to the electric dipole moments if the 
first two generations of sfermions are heavier than 0(1 TeV) Although somewhat tuned, these 

scenarios have sparked quite some interest, and need to be probed directly through collider searches of CP 
violation in the Higgs sector. 
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soft-breaking masses for the third generation squarks; M~ is the complex gluino mass, and 
ca-i c m are real numbers. (The mass parameters entering the slepton sector and the mass of 
the two weak gauginos are irrelevant for our discussion.) The remaining parameter needed 
to specify the Higgs sector is chosen here to be the pole mass of the charged Higgs boson, 
m H ± Strictly speaking the spectrum of CPX scenarios has ca = c M = 1. Here, we extract 
informations also for cases with ca, < 1. We take Msusy = 0.5 TeV; vary tan/3, the ratio 
of the v.e.v.'s and v u acquired by H% and H® at the minimum of the Higgs potential, and 
vary also . As for the phases [Ej present in these scenarios, we take as free parameters 
the arguments of the products A t /j,, A b fi, which we assume to be equal, and of the product 
Mgfi, respectively = arg(A t //) = aig(A b fi), and <£> 9A1 = arg(M^ / u). When relevant, we 
discuss also the four CP-conserving cases obtained for ^a^, ® gf i = 0°, 180°. Our numerical 
analyses make use of the recently-developed program CPsuperH [To] . 

We start observing that, once the threshold corrections to the 6-quark mass are included, 
the in-general-complex 6-quark Yukawa coupling is 

V2m b 1 

h = n , 2 

v cosp Kb 

where v 2 = v \ + t> 2 , with v ~ 254 GeV, and v d = v cos/3, v u — v sin/3. The factor R^. 

Rb = 1 + Kb tan (3 (3) 

collects in Kb the finite corrections § to the 6-quark mass. In turn, Kb can be split as 

K b = t g + e H , (4) 

where e g and e# are, respectively, the contribution coming from the sbottom-gluino exchange 
diagram and from the stop-Higgsino diagram. Their explicit form is 
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e * = 3# M ^* 7 K><> l M ?D > e " = 4V/K,<, W) • (5) 

The one-loop function I(a,b,c) is well known and can be found, for example, in [Hj. The 
left-right mixing elements in the matrices for the sbottom- and stop-mass squared used to 
obtain these corrections, are, in our convention, 

-^K( A b v d~ VV U ) , h* {A*v u - fiv d ) . (6) 

For the spectrum of Eq. (JTJ) with c A = c M = 1 and M SUS y = 0.5 TeV, it is |/z| 2 , |M~| 2 > M| USY 
and, as a consequence, the sbottom-gluino corrections are, in absolute value, much larger 
than those obtained from stop-Higgsino exchange. For a s and h t at the scale MgusY; i-e. 



$We remind that the issue of electroweak-symmetry breaking, possibly radiatively induced as in the 
constrained Minimal Supersymmetric Standard Model is not addressed in these scenarios. 

§ There are additional corrections Shi,/ hi, to Rb which are not enhanced by tan/?. We neglect them in our 
discussion but their effect is fully included in the numerical analysis jlfij . 
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a s ~ 0.1 and h t ~ 1, |e s | has a value ~ 0.05 and it is about one order of magnitude larger 
than |e#| (|e 9 | is larger than |e#| by the factor ~ 7r |//| 2 /|A t M~|) In this case, the threshold 
corrections to the 6-quark Yukawa coupling are strongly affected by the phase $ gM , whereas 
the dependence on is weak. We remind that, at the one-loop level, the corrections 
inducing the CP-mixing in the Higgs sector are sensitive only to $,4^. affects this 

mixing at the two- loop level through the one- loop corrections to the b- and t-quark masses |4, . 
Notice that the hierarchy \e g \ ^> \eh\ holds in general, except in the somewhat unlikely cases 
\M~\ <C | At | or | /x 2 1 , M| USY <C |A t M~|. Possible variations of the renormalization scale of ot s 
and ht do not affect this statement, neither have any substantial impact on the numerical 
results for the production cross sections of the various states H^. 

The effect of the radiative corrections becomes most significant for $ 5M = = 180°, 
when dramatic constraints on tan /3 can be obtained. Indeed, to prevent the lighter sbottom 
squark b\ from becoming tachyonic, the 6-quark Yukawa coupling is constrained as 



V \fi\ 



l + O 



\Ah 



\fi\ tan/3 



(7) 



where we have taken v u ~ v and neglected the m b - and D-term contributions to the diagonal 
elements of the matrix for the sbottom-mass squared. Since tan/3) is in general 

negligible, the constraint on h b can be easily recast into a constraint on tan j3. That is, the 
region of tan/3: 

' tan/3 < — (8) 



< 



m b \fi\ m b \jj\ 



is not allowed when 



/If 2 11 M 2 



mb\n 



\ Kb \ > TJ2 ' K b = -\n b \. (9) 

M SUSY 

These conditions are easily satisfied by the spectrum in Eq. (JTJ), for any value of ca and c M . 

The constraint of a nontachyonic b\ excludes the tan/3 region 12 < tan/3 < 30 for the 
scenario in Eq. (0) with M SU sy = 0.5 TeV, c A = = 1, and $ 9Ai = $ Afl = 180°. The 
value m b = 3 GeV was used here and will be used also for all other numerical evaluations 
throughout this paper. In such a case, the coupling h b has at most the value 0.7 in the 
allowed region, as shown by the left frame of Fig. By imposing that the lightest-Higgs 
boson mass exceeds, say, 115 GeV the excluded region gets extended up to tan/3 ~ 50. 
Notice that the limiting case tan (3 — > oo, i.e. Vd — > 0, is not a problematic one, when it comes 
to obtain acceptable values of mr and m Hl . The h b coupling is, in this case, h b ~ y/2m b /vn b 

"Such large values of \e g \ and \ch\ are specific to the scenarios considered here and are compatible with 
those found in Ref. |18| . 

"This value is reminiscent of the lower bound on the lightest-Higgs boson from LEP2 This was, 

however, obtained in the CP-conserving case and in the limit of a heavy CP-odd Higgs boson. Moreover, it 
is known JU] that in CP-violating scenarios the lower bound that can be deduced from LEP data may be 
considerably weaker than that reported in Ref. |19j . In our case case, the choice of this value is motivated 
by the fact that it allows simple comparisons with existing results in the literature. Particularly important 
will be later on the comparison with results presented in Ref. |2(J|. where the same choice is made. 
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Figure 1: Absolute value of the Yukawa coupling of the b quark, vs. tan/? (left frame) for the 



scenario 



of Eq. QJ with ca 



1, Msusy = 0.5 TeV, and <$> A u = $ 



ni-i 



180°. For the same 

scenario, but fixed values of^Afi, the regions of the plane (Qg^, tan/3) in which the mass of the b\ 
squark is vanishing are shown on the right frame. They are enclosed by a dashed, solid, dot-dashed 
line corresponding to <&Afi = 0°, 90°, 180°. 



and the mass of the b quark is generated at the quantum level J7j. For smaller values of 
c M , for example such that \fi\ = M SU sy, the excluded region of tan/3 at intermediate/large 
values is further enlarged by the fact that \hb\ grows rapidly above perturbative levels. 

This tan/3 exclusion occurs also for nontrivial values of $ 9M and ^a^-, as it can be seen 
in the right frame of Fig. ^ in it, we show explicitly the regions of the plane $ 9At -tan/3 in 
which the mass of the b\ squark is negative. They are enclosed by a dashed, solid, dot-dashed 
line, corresponding to — 0°, 90°, 180°, respectively. The scenario chosen for this frame 
is, again, that of Eq. (JTJ) with M SU sy = 0.5 TeV and ca = c M = 1. Regions of tan/3 to be 
excluded are found for $ 9A1 ~ 180° ± 30°. This pattern has to be compared with the very 
large tan /3 exclusion typical of constrained supersymmetric models with = $ 9At = 0°, in 
which the threshold corrections to mb are not taken into account [2T] . 

We turn now to consider the b-b-Hi couplings. Once the threshold corrections to the 
6-quark mass are included, the effective Lagrangian for the interaction of the neutral Higgs 
boson to b quarks can be written as 

C = - r ^b(gl + ig r lb )b^ (10) 

with = (</>i, 02, «), where 0i and 02 are the CP-even parts of the neutral components of 
the two doublets: 

H° d = -^(^ + 0i + m x ) , H° u = -= (v u + 02 + za 2 ) , (11) 

and a is a combination of the two CP-odd components a\ and 02, a = — ai sin/3 + 02 cos /3. 
The orthogonal combination, G° = a\ cos /3 + 02 sin (3, is the Goldstone mode. The couplings 
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Figure 2: Couplings g^ P vs. <3? SAt , for the spectrum in Eq. with ca = = 1, Msusy = 0.5 TeV, 
and tan/3 = 10. The solid lines are for ^a^ = 0°, the dashed ones for §Afi = 180°. The horizontal 
lines indicate the values of the uncorrected couplings. 
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S,P 



are [THj 



1 Stefi-V a* = ^^ Kb 



cos P \RbJ 01 cos P \R, 



■ s -— *>(£). «£ = -^ (12) 



:>2 



cos f3 \RbJ ' '"' cos P \Rb< 

( tan (3 — Kb 



g s a = (tan 2 /3 + l)3m(^), ^ 



\ Rb 



and for values of tan/? such that |/«jt,| tan/5 ~ 1 (with one of the two, or both possibilities: 
$le(Kb) tan/3 ~ 1, tan /3 ~ 1) reduce to 

c tan 3 r <y . , v 01 p tan /3 r . . - 

% = r^ja i 1 + yte(Kb) tan 3\ , = — — [3ta(K 6 ) tan/3] , 

9te = TETn [^e(K 6 )tan/3+ |k 6 | 2 tan 2 /?] , = —rig [Sm(/e&) tan/3] , (13) 

\tib\ L J |-n*r 

= TTTif [9to(«6)tan/3] , #f = [1 + 3?e(ft 6 ) tan/3] . 

Before proceeding further, we list in the following some of the interesting features of these 
couplings. We illustrate them numerically for a specific CPX spectrum, i.e. with ca = c M = 1, 
-^susy = 0.5 TeV, and tan /3 = 10, which is a value not plagued by the problem of a tachyonic 
b\ squark discussed earlier. For visual clarity we also illustrate some of these features in Fig.EJ 

• If no threshold corrections are included, the only nonvanishing couplings are = 1/ cos /3 
and g£ = — tan/3. The inclusion of these corrections affects these two couplings mainly 
through the factor 9?e(l/i4), which is a suppression or an enhancement factor, depending 
on the value of and varies between ~ l/(l + |e 9 | tan/3) and ~ 1/(1 — |e fl | tan/3), obtained 
for $ 9Ai = 0° and $ 5Ai = 180°, respectively. We have used here the fact that \e g \ ^> |e^|, 
which, as previously discussed, is rather generic in our scenarios. Notice also that the factor 
9fte(l/.Rfc) is larger than 1 for cos $ gA1 < — |e ff | tan/3. For tan/3 = 10 and the spectrum specified 
above, this happens when 135° < $ 9A1 < 225°, and and \g£\ reach the maximum of 15 
or ~ 20, depending on the value of <& AfJ , (i-e. 0° or 180°), at $ SAt = 180°. 

• The inclusion of threshold corrections is responsible for the appearance of the other four 
couplings: g^ 2 and g^ 2 , the smallish ones, i.e. without the overall factor of tan/3 that " s 



g£ have (when |k 6 | tan/3 ~ 1), g>£ and , the large ones, with the tan/3 factor. 

• If no CP phases are present, only gf 2 , and g£ are nonvanishing and, again, the effect 
of the threshold corrections is more significant for Q gfl = 180° than for $ 9M = 0°. 

• As the phase varies, the moduli of the couplings g? v g^ 2 , and gf reach their extremal 
values at cos$ 9M ~ — 2\e g \ tan/3/(l + |e g | 2 tan 2 /3). In the specific case considered here, this 
is 5 or 8 for \g?\ and |gf| at Q gfJi ~ 150° and 210°, depending on the value of In 
contrast, \g% 2 \, as well as \g% 2 \, never exceed 1.2. 

The states <f>\, 02, and a are not yet the neutral Higgs boson mass eigenstates. Their real 
and symmetric 3x3 matrix M. 2 H has nonvanishing entries that mix the two states 4>\ and 4>2 
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as well as nonvanishing CP- violating entries that mix a with <j>\ and 02 and are proportional 
to &Afi- The diagonalization of this matrix through an orthogonal matrix O, 



T M 2 H = diag 



m H 1 -> m H 2 -> " H 3 J > 



(14) 



yields the three eigenstates H±, H2, H3, ordered for increasing value of their masses. In the 
limit of vanishing Hi is h, H2 and H% are H and A, or vice versa, A and H, depending on 
the values of rriA and tuh- After this rotation, the effective Lagrangian of Eq. fllOl) becomes 



where and are 



Oai gi 



p — n p 

9 Hi — U ai 9a J 



(15) 
(16) 



with the index a running over (0 1; </> 2 , a) and i over (1, 2, 3). 

We are now in a position to evaluate the production cross sections of the neutral Higgs 
bosons Hi via 6-quark fusion at hadron colliders. They can be expressed as: 



<j(hadihad2 — > bb — > Hi) = cx(bb — > Hi 



dx 



— b hadl (x, Q) 6 had2 ( — , Q) + (6 <-»• 6) 



(17) 



where &hadi(^> Q) and &had,(^, Q) are the 6- and 6-quark distribution functions in the hadron 
hadj, Q is the factorization scale, and the Drell-Yan variable n = m H ./s, with s the 
invariant hadron-collider energy squared. Finally, the partonic cross section is 



a(bb -> fli) 



The sum of the couplings squared 



v 2 6m H . 



(9nf + G&) 



P \2 



reduces, in the limit in which no threshold corrections are included, to 



{9 S H f + (9h, 



P \2 



no— thresh— corr 



P\2 



(19) 



(20) 



with g^ = 1/ cos/3 and g£ = —tan/5. For i — 1, in the limit of large m#±, O^i — > cos/5 
and O a i — > 0, leading to the usual Standard Model coupling of the lightest neutral Higgs 
boson to b quarks, in both, CP-conserving and CP-violating scenarios. That is, the overall 
tan/? dependence of the coupling <?? is killed by the mixing element 0^\. 

Through the inclusion of the threshold corrections to the b quark in CP-conserving sce- 
narios, the element S e ^ s switched on. Being one of the small couplings, i.e. one without 
an overall tan/? factor when tan/3 is large, g^ 2 is not expected to produce great numerical 
deviations in the values of the cross sections at large tan/3. Deviations of 0(1) can only be 
induced by the factors l/|i?b| 4 , which enhance or suppress the uncorrected cross sections, 
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depending on the sign of the threshold corrections. It is however in CP-violating scenarios, 
when two couplings with an overall tan (3 dependence are switched on, gf and gf , that 
the pattern of Higgs production through 6-quark fusion can become very different. This, 
of course, if the projection of Higgs boson current-eigenstates to mass-eigenstates is not 
particularly destructive. 

To analyze the behaviour of the sum in Eq. ([EH) , we use the approximation: 

[& + (g p H f 

In this, we have parametrized the sum [{gfj-) 2 + (fi^) 2 ] in terms of the only elements O^i 
and O a i that appear in Eq. (|2Uj) . i.e. in the case in which no threshold corrections to the 
6-quark mass are included. This approximation relies on the properties of the couplings g^ P 
discussed above, i.e. on the fact that g^ ~ —g p , g P x ~ <?f , and that g^ 2 and g P 2 can be 
neglected. In the scenarios considered here, this approximation turns out to be valid for 
most of the relevant values of the phases <3> 9/i and ^a^- 

To strengthen our point, we show explicitly the mixing elements for the scenarios 
already considered in Fig- El We show them in Fig. Inversus $a^, but with <3> 5M fixed at 180°. 
As already remarked, the elements have a rather weak dependence on $ gM , coming from 
the two-loop corrections to the Higgs potential. The charged Higgs boson mass is solved in 
order to have m Hl = 115 GeV, for all values of an d and it is therefore different for 
different values of these phases. As a consequence, also m H2 and m^ 3 are varying quantities. 
Notice that, for 5^ « 100°, H% is predominantly the CP-odd a boson, whereas H 2 and 
H3 are mainly 02 and <pi, respectively. For these values of Qa^, the charged Higgs boson is 
relatively light, i.e. it has a mass < 200 GeV. 

A comparison of Figs. 121 and El indicates that the approximation in Eq. (J2TJ) is probably 
not adequate when [0\ xi + 2 J <C 1, since the neglect of the elements g^ 2 and g? cannot 
be justified. Therefore, we keep the exact expression in Eq. (|19|) for the calculation of 
the production cross sections, but we use this approximation in order to clarify the tan/3 
dependence of the sum 

[(9hJ 2 + (9 Hi) 2 ]- % substituting in Eq. (0} the expression of the 

S P 

couplings g± at large tan/3, we obtain: 

+ - (°5x* + °y^r- ( 22 ) 

That is to say, the terms in tan 4 /3 cancel identically. It remains only a factor of tan 2 f3 sup- 
pressed or enhanced by l/|i?;,| 2 , in which all dependence on <3> 9M gets concentrated. Whether 
this factor of tan 2 f3 is inherited unharmed by the cross sections for all Higgs bosons Hi 
depends on the matrix elements O ai and 0^. For the values of $^ for which the charged 
Higgs boson is not very heavy, and the three Higgs boson states are significantly mixed, we 
expect this to be the case, in general for all three states, and in particularly also for the 
lightest one. That is, for $^ w 100°, we expect a large enhancement of the production 
cross section for Hi, but presumably a suppression of that for H 2 , as an inspection of Fig. El 
suggests. 
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(21) 
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Figure 3: JTie sums 0?^ + us. ^a^, for the spectrum of Eq. £7J} with ca = c M = 1, tan /3 = 10, 
as in Fig.ffl and mjf a = 115 GeF, <& gtx = 180°. T/ie dashed lines show O ai ; the dotted ones, O 2 ^. 



We show the partonic cross sections in Fig. as functions of $A/x, for <E> ffM = 0° (dashed 
lines) and 180° (solid lines). The dotted lines indicate the cross sections without threshold 
corrections to the 6-quark mass. The cross sections for all neutral Higgs bosons are shown for 
tan p = 10, that for Hi also for tan /3 — 5. As mentioned, the two upper frames show cross 
sections for constant uihu i-e. = 115 GeV, which cannot actually be realized for all 

values of <3>^ M in the case of tan /3 = 5 (m Hl tends to be lighter for < 50° and ^a^ > 130°). 
In contrast, in the two lower frames, m# 2 and mn 3 are different for different values of 
mH 2 (. m H 3 ) reaches the maximum of ~ 240(250) GeV at = 0° and $ 9A1 = 180°, the 
minimum of ~ 120(150) GeV for 3^ ~ 90°, nearly independently of <3> 9M . All cross sections 
show the typical pattern already observed for 0?^ + 2 ai in Fig. 01 with modulations in the 
cases i = 2, 3, due to varying values of the corresponding Higgs masses. Notice the increase 
of almost two orders of magnitude in the cross section for Hi, at $^ ~ 100°, when the 
value of tan/3 is only doubled. Indeed, for i — 1, the sum 0?^ + 2 ai is larger for tan/3 = 10 
than for tan/3 = 5, when larger values of are needed to ensure that = 115 GeV. 
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Figure 4: Partonic cross sections for the b- quark fusion production of Hi, H2, and H3 vs. &A[i, 
for $> giM = 0° (dashed lines) and 180° (solid lines). The supersymmetric spectrum is that used also 
for Fig. In dotted lines are also the cross sections with no nib- corrections. 

Moreover, at tan/3 = 5 the eigenstate Hi has still a large CP-even component. 

Overall, the mj corrections increase the cross sections a(bb — > Hi) with respect to the 
uncorrected ones for = 180°, decrease them for = 0°. The impact of these correc- 
tions, is in general always large, producing enhancements up to one order of magnitude and 
suppressions down to —60%. Notice that the corrections remain surprisingly large for H 2 and 
H 3 , also in the CP-conserving cases $ 3M = 180°, $ Afl = 0° (+120%); $ 9M = 180°, $ Am = 180° 
(+300%); as well as $ 9M = 0°, = 0°, 180° (-50%). They are more modest for Hi at 
$ 9M = 0°, 180°, = 0° (i.e. ~ ±20%), but still +70% and -50% at = 0°, 180° and 
= 180°. 

After convoluting the parton distribution functions, we obtain the hadronic cross sections 
for the Tevatron (y/s = 1.96 TeV) and the LHC (\/s = 14TeV). These are shown in Fig. El 
vs. §A[i, for tan/3 = 10 and two values of $ 9A1 = 0° (dashed lines) and $ SM = 180° (solid 
lines). We have used the leading-order CTEQ6L [22] parton distribution functions and chosen 
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the factorization scale Q = m^/4:. This has been suggested in most of the papers in Ref. [23] 
as the scale that minimizes the next-to-leading-order QCD corrections to these cross sections 
when no threshold corrections to are kept into account. Although this should be explicitly 
checked, we believe that the inclusion of these corrections, which amounts to substituting the 
tree-level Yukawa couplings with effective ones, should not affect substantially this result. 
We notice also that these supersymmetric threshold corrections capture the main part of all 
supersymmetric corrections to the production cross sections of neutral Higgs bosons through 
6-quark fusion. Other corrections, with a nontrivial dependence on the momenta of the Hi 
bosons are of decoupling nature, and therefore subleading. (See the explicit check for the 
related processes gb — > bHi in the CP-conserving case of Ref. [24J.) 

As already observed at the partonic level, all three production cross section can deviate 
considerably from those obtained in CP-conserving scenarios and the impact of the super- 
symmetric threshold corrections to m;, is very important. The production cross section for 
the lightest Higgs boson through the 6-quark fusion is, in general, comparable to the pro- 
duction cross section via gluon fusion. For some values of $ 9A1 and Qa^, it can even be 
larger. For example, for $ 9A1 = ~ 100°, a production cross section via gluon fusion of 
< 30 pb at the LHC, which can be easily evinced from Ref. [2U| for the same set of super- 
symmetric masses used here, in particular the same value of rriH 1 , is indeed smaller than 
the corresponding cross sections in Fig. (The analysis of Ref. (201 does not include the 
supersymmetric threshold corrections to m&. Their impact, however, is expected to be less 
dramatic in the case of the gluon- fusion production.) Even for different values of $ 9A1 and 
$Afi, the 6-quark fusion production mechanism cannot be easily dismissed as subleading. 
No comparison between the two production mechanisms for the two heavier Higgs bosons 
is possible at the moment, since in the existing studies [25J of H2 and H3 production via 
gluon fusion different supersymmetric spectra than those considered here are analyzed. For 
all three mass eigenstates Hi, however, the cross sections corresponding to the two produc- 
tion mechanisms have rather different dependences on the supersymmetric spectrum. The 
results shown here for the 6-quark fusion cross sections, for example, remain unchanged for 
values of M SU sy different from our representative value of 0.5 TeV, provided the relative size 
of the different supersymmetric parameters is not changed. The same is not true for the 
gluon-fusion cross sections, which are sensitive to the absolute value of Msusy- 

More investigations, theoretical and experimental, are needed to unravel all implications 
of the enhanced cross sections presented here. A detail comparison of the yield of all neutral 
Higgs bosons of supersymmetric scenarios through 6-quark fusion, as opposite to the yield 
through gluon fusion and through Higgs strahlung (which is relevant for the Tevatron), will 
have to be performed. Moreover, all subsequent decay modes for these Higgs bosons need 
to be studied to allow unambiguous interpretations of possible Higgs bosons signals, which 
hopefully will be detected at the Tevatron and/or the LHC. All future analyses need to be 
carried out at the same level of precision, i.e. incorporating the threshold corrections to m& 
that, as proven in this paper, turn out to be very important in these scenarios, not only for 
large but also for intermediate values of tan (5. 

To summarize, we have studied the effects of threshold corrections to the 6-quark mass 
on the production of neutral Higgs bosons via 6-quark fusion, bb — > Hi, in supersymmetric 
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Figure 5: Cross sections for the production of Hi, Hi, and H3 vs. <&An, for <& gil = 180° (solid 
lines) and 0° (dashed lines), at the LHC with y/s = lATeV (two upper lines) and at the Tevatron 
with y/s = 1.96 TeV (two lower lines), hadihacb = pp for the LHC, pp for the Tevatron. 



scenarios with large CP-violation in the Higgs sector. This is assumed to be induced by 
explicit phases in super symmetric and supersymmetry-breaking parameters. For these par- 
ticular scenarios, we have found that: 

i) large phases of the combination Mgfi, i.e. $ 9M = 180° ± 30°, for all values of can 
drive the mass squared of the lightest b squark to negative values and, depending on the 
supersymmetric spectrum, the 6-quark Yukawa coupling to nonperturbative ones; 

ii) large deviations in the behaviour of the 6-quark fusion cross sections, with respect to the 
same cross sections in CP-conserving scenarios, are obtained for ~ 100° when the Higgs 
mixing is maximal; 

Hi) the supersymmetric corrections to the 6-quark fusion production cross sections are in 
general very large, even for the four CP-conserving scenarios obtained by fixing the values of 
$ 9M and &A/i at 0° or 180°. Among these four scenarios, those with <3> ffM = & A/X = 180° give 
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the largest cross sections cr(hadihad2 —>■ bb — > Hi) both at the Tevatron and at the LHC. 

We conclude that in CP-violating supersymmetric scenarios, the production of neutral Higgs 
bosons through 6-quark fusion cannot be neglected, in general, with respect to the produc- 
tion through gluon fusion. Dedicated theoretical and experimental studies should be carried 
out in order to allow a reliable detection of these mixed state neutral Higgs bosons, or to 
constrain the CP- violating scenarios that predict them. 
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